Successive transitions and intermediate chiral phase 
in a superfuilid 3 He film 
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Superfluidity ordering of thin 3 He films is studied by Monte Carlo simulations based on a two- 
dimensional lattice spin model with Z2 x U{1) x 5*0(3) symmetry. Successive phase transitions 
with an intermediate 'chiral' phase, in which the Z-vector aligns keeping the phase of the condensate 
disordered, is found. Possible experiments to detect the successive transitions are discussed. 

67.70.+n, 67.57.Lm 



It is now well established that the superfluidity transi- 
tion of 4 Hc films belongs to the Kosterlitz-Thouless (KT) 
universality class governing a variety of two-dimensional 
(2D) systems with U(\) symmetry [1]. By contrast, 
the nature of the 2D phase transition of its isotope, 
i.e., the superfluidity transition of 3 Hc films of thick- 
ness of order the zero-temperature coherence length, 
0.01/im IS Co ~ 0.05/im, has not been well understood. 
Although the complex structure in the order parame- 
ter of 3 He [2] makes its transition behavior an intrigu- 
ing problem, the experimental observation of the possible 
nontrivial critical behavior of bulk 3 He is difficult, since, 
due to its long coherence length, the critical regime is 
limited to a narrow temperature range close to T c , i.e., 
t =\ (T — T c )/T c \< t* k 10~ 6 [3]. By contrast, in the 
case of thin films, enhanced effects of fluctuations due to 
the low dimensionality give rise to much wider critical 
regime, as wide as t* w 1CP 1 ~ 10~ 2 [4]. Hence, there 
appears to be a good possibility to experimentally ob- 
serve the novel transition behavior originated from the 
complex order parameter. 

There have been many experimental works on the su- 
perfluidity ordering of 3 He in restricted geometries in- 
cluding cylindrical pores [5,6] and thin flat films [7,8]. 
Indeed, several experiments reported the observation of 
novel ordering behaviors, even a signature of new phase 
[6,7]. Experimentally, one challenging problem is to make 
measurements on a sample which can be regarded as a 
good 2D system, i.e., a thin flat film of uniform thickness 
of order 0.1/im or less on sufficiently specular boundaries. 
For such an ideal 2D system, theoretical calculation sug- 
gests that the stable superfluidity state is that of the A 
phase down to zero pressure [9] , which has also been con- 
firmed experimentally [8]. The order parameter can be 
described as a triad of three orthogonal unit vectors m, 
n and 1 = m x n, where the 1-vector is kept perpendicular 
to the surface. 

Some time ago, Stein and Cross argued that such su- 
perfluidity state with the A-phase structure possessed the 
novel twofold Ising-like degeneracy, according as the 1- 
vector was either parallel or antiparallel with the surface 
normal [10]. This Z2 symmetry coexists with the stan- 
dard U(l) gauge symmetry associated with the phase of 
the condensate. The latter symmetry corresponds to a 



continuous rotation of the m- and n-vectors around the 
1-vector. Stein and Cross then argued that this novel Z 2 
symmetry might lead to an Ising-like transition which 
might "compete" with the standard KT transition as- 
sociated with the U(l) gauge symmetry, although the 
detailed nature of the transition was not specified. 

Since then, 2D phase transitions associated with the 
Z2 x U(l) symmetry have been studied in a different 
area, i.e., in the context of frustrated 2D XY systems 
such as the triangular-lattice XY antiferromagnet [11] 
or the Josephson-junction array in a magnetic field [12]. 
In these problems, the Z2 Ising-like degree of freedom 
has been called 'chirality' [13]. While these studies have 
established the occurrence of a phase transition with a 
sharp specific-heat anomaly, there still remains some con- 
troversy as regards whether the Z(2) and U(l) degrees of 
freedom order simultaneously, or at two close but distinct 
temperatures. Meanwhile, some differences exist between 
these frustrated XFsystems and 3 He films. Particularly 
important point neglected in the previous analysis [10] 
may be the fact that 3 He has an additional internal sym- 
metry, SO (3), associated with the nuclear spin degree of 
freedom of the condensate. 

In the present Letter, I perform a numerical study 
of the superfluidity ordering of 3 He films in order to 
clarify possible thermodynamic phases and the transi- 
tion behavior between them. For that purpose, I intro- 
duce a 2D lattice spin (pseudospin) model which pos- 
sesses the expected full symmetry of the order parameter, 
Z2XU (1) x SO(3), and perform Monte Carlo simulations. 
It is found that the model exhibits two successive transi- 
tions with a 'chiral' intermediate phase. 

The order parameter describing the 2D superfluidity 
of 3 He may be described by a 3 x 2 tensor variable Anj, 



A^j = Ad^mj + irij), 
1 = 1 n 1= 1, m n = 0, 



(1) 



where fi = x,y, z refers to the spin component, j — x,y 
refers to the real-space coordinate (the film surface is 
taken to be the xy-plane) , d is a three-component unit 
vector in spin space representing the spin state of the 
condensate, while m and n are mutually orthogonal two- 
component unit vectors in real space representing the or- 
bital state of the condensate. The Ising-like variable, or 
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the chirality, is defined by r = sign(7 z ) = m x n y — m y n x , 
where t = ±1 represents the 1- vector pointing either up 
or down along the surface normal. 

In deriving an appropriate model Hamiltonian, I start 
with the standard Ginzburg-Landau free energy in the 
London limit, 

H GL « K{\ VA M | 2 +S | divA M | 2 ), (2) 

where (A li )i=A li i. Although there generally exists a 
spatially-anisotropic gradient term (the second term of 
Eq.(2)), rcnormalization-group e = 4 — d expansion anal- 
ysis showed that such spatial anisotropy was irrelevant at 
the superfluidity transition, i.e., S — ► upon renormal- 
ization [3] . Hence, for simplicity, I drop the second term 
here. Neglecting the weak dipolar interaction of order 
10~ 1 /iK, restricting the space to 2D, parametrizing the 
m-vector as m = (cos 8, — rsin#) where < 6 < 2n is 
the phase angle of the condensate, and discretizing the 
continuum into the lattice, one obtains 

U = - J J2 (1 + TiTj) COs(0i - 6j)di ■ dj, (3) 

<ij> 

where J > is a coupling constant and the summation 
is taken over all nearest-neighbor pairs on the square lat- 
tice. The Hamiltonian (3) can be viewed as a coupled 
Ising-Xy-Hcisenberg model with Z 2 x U(l) x 50(3) sym- 
metry, in which Tj = ±1 is an Ising variable [Z 2 ], 0i be- 
ing an angle variable of the Xy-pseudospin pi = 
(cos#i, sin#,), and is a Heisenberg variable [5*0(3)]. 
Note that the Hamiltonian (3) also possesses the local 
symmetries under the transformations; (i) pi — > npi; 
(ii) di —> ndi] (hi) pi -> ±p; and d { -> ±d;. 

The ordering of the 1-vector, or of the chirality, can 
be probed via the Ising-magnetization, r = n (N 

is the total number of lattice sites), by calculating the 
average chirality and the associated Binder ratio, 

^<- 2 > V2 , 9r = \{*-^). (4) 

Since the correlation functions associated with the phase 
variable 8, or the Xy-pscudospin p, are not invariant 
under the above local transformations and vanish triv- 
ially, one needs to calculate the "local-gauge-invariant" 
quantity to study the phase ordering. An appropri- 
ate quantity is a second-rank symmetric traceless ten- 
sor P illv = -^(p^Piv - \5^ v ), or equivalently, a new 
Xy-pseudospin variable with the doubled phase angle, 
p- = (cos2#j,sin2#j). Via the Xy-magnetization, p' = 
lv Si Pi' ^ e associated Binder ratio is defined by 



Likewise, the ordering of the d-vector is probed via 
= jjJ2iDifiu, where a second-rank symmetric 



traceless tensor is defined by Z?i M „ = -^(di^d^ — ^S^). 
The corresponding Binder ratio is given by 

^ = ^-5^2), ^ 2 -E^- (6) 

flL* 

Monte Carlo simulations are performed based on the 
standard single-spin-flip Metropolis method. The lattice 
contains N = L 2 sites with L = 20,30,40 and 60 with 
periodic boundary conditions. Typically, I generate total 
of 10 5 Monte Carlo steps per spin (MCS) at each tem- 
perature, and (2.5 ~ 10) x 10 5 MCS in the transition 
region. 

The calculated specific heat is shown in Fig. 1 . A sharp 
peak is observed at 77 J = 0.754 ± 0.002. As is shown 
in the inset, the peak grows slowly with increasing L, in- 
dicative of a continuous transition characterized by the 
exponent a ~ 0. 
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FIG.l The temperature and size dependence of the spe- 
cific heat per site. The inset represents a magnified view 
of the transition region. 

The calculated Binder ratios (4)- (6) are shown in 
Figs.2a-c, respectively. As can be seen from Fig. 2a, g T 
for different L cross at a temperature which coincides 
with the specific-heat-peak temperature. This indicates 
the occurrence of a continuous transition at T c \/J = 
0.753 ± 0.001 into the phase with a finite chiral long- 
range order, f > 0, accompanied by a sharp specific-heat 
anomaly. The behavior of f (not shown here) turns out 
to be fully consistent with this. By contrast, as can be 
seen from Figs. 2b and c, gg and gd for various L donot 
cross or merge at T c \ , constantly decreasing with increas- 
ing L around T c \, a behavior characteristic of a disorder 
phase. Hence, the state just below T c \ is a pure chiral 
state with only a chiral long-rage order while the phase 
and the spin are kept disordered. 

At a lower temperature T c2 /J = 0.692 ± 0.04, gg for 
various L merge, and below T < T c2 , continue to stay on 
a common curve: See Fig. 2b. Such a behavior is 
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FIG. 2 The temperature and size dependence of the 
Binder ratios in the transition region; (a) of the r- variable 
(chirality); (b) of the p'-variable (phase); (c) of the d- 
variable (spin). Note the difference in the temperature 
scale in each figure. The insets show log-log plots of the 
L dependence of the ordering susceptibilities, N < t 2 > 
(a), N < p' 2 > (b) and N < D 2 > (c). In the ordered 
phase with a nonzero order parameter, the data for larger 
L should lie on a straight line with a slope equal to two, 



while in the disordered phase the data should exhibit a 
downward curvature staying finite even in the L — > oo 
limit. At criticality and in the KT-like phase, the data 
should lie on a straight line with a slope 2 — 77, which is 
less than two. 

characteristic of the KT-like transition in which the low- 
temperature phase is a critical phase with algebraically- 
decaying correlations. Thus, the quasi-long-range order 
of the £7(1) phase is realized below T c2 . The absence 
of any appreciable anomaly in the specific heat around 
T c2 is also consistent with such KT-like transition. By 
contrast, monotonically decreases with increasing L 
around T c2 , suggesting that the d- vector remains disor- 
dered even at T < T c2 . This observation is consistent 
with the fact that the d-vector is a Heisenberg variable 
carrying the SO(3) symmetry, and that the 2D Heisen- 
berg model orders only at T = 0. I have also checked 
that the behaviors of the order parameters corroborate 
the above conclusion. As an example, the L dependence 
of the ordering susceptibilities, N < r 2 >, N < p' 2 > 
and N < D 2 >, are shown on log- log plots in the insets 
of Fig.2. 

Thus, the present model exhibits two successive tran- 
sitions at T = T cl and T c2 , where T cl ~ 0.753 J > T c2 ~ 
0.692J. The intermediate phase is a pure chiral state 
in which only the Z 2 chirality, representing the direction 
of the 1- vector, aligns. In the low-temperature phase, in 
addition to the chiral long-range order, the U(l) phase 
of the condensate exhibits the KT-type algebraic order 
leading to a true superfluidity. Meanwhile, the d-vector 
representing the spin state of the condensate remains 
disordered at any finite temperature. Although in real 
films the dipolar interaction neglected here gives rise to 
additional interaction, it is several orders of magnitude 
weaker than J [2] and the ordering of the d-vector, if any, 
should occur at a temperature much lower than T c i and 
T c2 . 

To identify the universality class of each transition, 
standard finite-size scaling analysis is made. The esti- 
mated chirality exponents at T = T c \, v = 1.0 ± 0.1 and 
(3/v = 0.13 ± 0.01, agree well with the values of the 2D 
Ising model, and I conclude that the upper transition 
belongs to the standard 2D Ising universality class. As 
regards the lower transition, the decay-exponent r/g as- 
sociated with the correlation function of the p'-vector 
at T = T c2 , < p'(r) • p'(0) >« r _w , is estimated 
from the log-log plot of N < p' 2 > versus L, yielding 
r]e = 0.45 ± 0.03: See the inset of Fig2b. The estimated 
rjg is considerably larger than the standard KT value, 
rj = 0.25, suggesting that the lower transition is not of 
the standard KT universality. It should be noticed, how- 
ever, that the spatial anisotropy neglected in deriving our 
model Hamiltonian (3) might possibly affect the value of 

Ve- 
in the present model, the intermediate phase is realized 
over an appreciable temperature range, in sharp contrast 
to the case of frustrated XY models where T cl and T c2 



3 



are identical or very close [11,12]. This difference is prob- 
ably caused by the SO (3) degree of freedom associated 
with the d-vector. Namely, although the d-vector itself 
does not order at any T > 0, its existence considerably 
affects the effective interaction between the chiral and 
phase variables. This may be seen by deriving an effec- 
tive Hamiltonian TL e ff written in terms of the r- and 
^-variables, by integrating out the d-variable. At suffi- 
ciently high temperatures, one can explicitly carry out 
this procedure to get 

H eff /k B T w - i K ^ 1 + cos(20j - 28 3 ) 

<ij> 

+K 2 T i T j }, (7) 

where one has K\ = K 2 . Indeed, a numerical study in 
Ref. [14] suggested that in the case of K i = K 2 the model 
exhibited two successive transitions with an intermediate 
chiral phase consistent with the present result, whereas 
the frustrated XY models corresponded to K 2 — where 
the model (7) was estimated to exhibit a single transition 
only. 

Finally, I discuss experimental implications of the ob- 
tained results. Main conclusion of the present work is 
that successive transitions and an intermediate phase are 
likely to be realized in thin films of 3 He. While the clear- 
est sign of the upper transition is a logarithmically di- 
vergent anomaly in the specific heat, it might practically 
be difficult to get the necessary experimental sensitiv- 
ity because specific-heat measurements usually require 
a considerable mass of sample. Some anomaly might 
be detectable at T = T c i, however, in the quantities 
like thermal-expansivity, optical properties, and sound- 
velocity or attenuation (particularly, the third sound 
[15]). By contrast, the NMR shift, which has widely been 
used in identifying the bulk superfluidity transition, is 
blind to the distinction between 1 and —1, and may not 
be a very good indicator of the transition. Rather, the 
shift should set in somewhat above T c \, by an amount 
of order t*, when the short-range order is developed and 
the 1-vector stays perpendicular to the surface. 

Although the intermediate phase is not a true super- 
fluid, it is expected to be a low-dissipation state distin- 
guishable from the high-dissipation normal state above 
T c i. Interestingly, the observation of a low-dissipation 
state was reported in the literature [5,6]. By contrast, 
the lower transition is characterized by the onset of a true 
superfluidity, and is detectable by the standard torsional- 
oscillator measurements. While the areal superfluid den- 
sity p s is expected to show a non-KT jump at T = T c2 , 



Ps/T c2 = 8m 2 kB/{^r)gh 2 ) (m is the mass of 3 He atom), 
its absolute value p s /T c2 ~ 17.6 x 10~ 12 g/(cm 2 mK) is 
small reflecting the low T c of 3 He, and experimental ob- 
servation of the jump itself might be difficult. 

Anyway, in the 2D superfluidity ordering of 3 He, new 
interesting physics different from that of helium four is 
expected. Further experimental studies are encouraged. 
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University of Tokyo. 



[I] J.M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 
(1973); D.R. Nelson and J.M. Kosterlitz, Phys. Rev. Lett. 
39, 1201 (1977). 

[2] A.J. Leggett, Rev. Mod. Phys. 47, 1159 (1977). 
[3] D.R. Jones, A. Love and M.A. Moore, J. Phys. C9, 743 
(1976). 

[4] One can roughly estimate t* based on the Ginzburg crite- 
rion, t* ~ [(27r^,) _d /AC] 2/(4 " d) , where d is the dimension- 
ality and AC is the specific-heat jump at T c [Modern Theory 
of Critical Phenomena, S.K. Ma, Benjamin, p. 94 (1976)]. If 
one puts d — 2 and substitutes typical values for the A-phase 
of 3 He, one gets t* ~ 0.04. 

[5] M.T. Manninen and J.P. Pekola, J. Low Temp. Phys. 52, 
497 (1983). 

[6] V. Kotsubo, K.D. Hahn and J.M. Parpia, Phys. Rev. Lett. 
58, 804 (1987). 

[7] J. Xu and C. Crooker, Phys. Rev. Lett. 65, 3003 (1990). 
[8] MR. Freeman and R.C. Richardson, Phys. Rev. B41 
(1990) 11011. 

[9] Y.-H. Li and T.-L. Ho, Phys. Rev. B38 (1988) 2362. 
[10] D.L. Stein and M.C. Cross, Phys. Rev. Lett. 65, 504 
(1979). 

[II] S. Miyashita and H. Shiba, J. Phys. Soc. Japan 53, 1145 
(1984); D.H. Lee, J. D. Joannopoulos, J.D. Negele and D.P. 
Landau, Phys. Rev. Lett. 52, 433 (1984). 

[12] S. Teitel and C. Jayaprakash, Phys. Rev. B27, 598 (1983); 
G. Ramirez-Santiago and J.V. Jose, Phys. Rev. B49, 9567 
(1994); P. Olsson, Phys. Rev. Lett. 75, 2758 (1995). 

[13] H. Kawamura, J. Phys. Condens. Matter 10, 4707 (1998). 

[14] J. Lee, E. Granato and J.M. Kosterlitz, Phys. Rev. B44, 
4819 (1991). 

[15] T. Ohmi, M. Nakahara, T. Tsuneto and T. Fujita, Prog. 
Theor. Phys. 68, 1433 (1982). 



4 



